Bose-condensed Bright Solitons under Transverse Confinement 
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We investigate the dynamics of Bose-condensed bright solitons made of alkali-metal atoms with 
negative scattering length and under harmonic confinement in the transverse direction. Contrary to 
the ID case, the 3D bright soliton exists only below a critical attractive interaction which depends 
on the extent of confinement. Such a behavior is also found in multi-soliton condensates with box 
boundary conditions. We obtain numerical and analytical estimates of the critical strength beyond 
which the solitons do not exist. By using an effective ID nonpolynomial nonlinear Schrodinger equa- 
tion (NPSE), which accurately takes into account the transverse dynamics of cigar-like condensates, 
we numerically simulate the dynamics of the "soliton train" reported in a recent experiment (Nature 
417 150 (2002)). Then, analyzing the macroscopic quantum tunneling of the bright soliton on a 
Gaussian barrier we find that its interference in the tunneling region is strongly suppressed with 
respect to non-solitonic case; moreover, the tunneling through a barrier breaks the solitonic nature 
of the matter wave. Finally, we show that the collapse of the soliton is induced by the scattering on 
the barrier or by the collision with another matter wave when the density reaches a critical value, 
for which we derive an accurate analytical formula. 
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I. INTRODUCTION 

The experimental achievement of Bose-Einstein con- 
densation with alkali-metal atoms at ultra-low tempera- 
tures [1,2] has opened the exciting possibility of studying 
topological configurations of the Bose-Einstein conden- 
sate (BEC), like solitary waves (solitons) with positive or 
negative scattering length a s [3]. Dark solitons (a s > 0) 
of Bose condensed atoms have been experimentally ob- 
served few years ago [4], while bright solitons (a s < 0) 
have been detected only very recently with Li using an 
optical red-detuned laser beam along the axial direction 
of the sample to impose a transverse (radial) confinement 
[5]. 

Bright solitons have been studied in one-dimensional 
(ID) models [6] and no investigation has been performed 
on the effects a finite transverse width which is always 
present in experiment. In this paper the dynamics of 
a Bose condensate is investigated by using the full 3D 
Gross-Pitaevskii equation [7] and an effective ID non- 
polynomial nonlinear Schrodinger equation (NPSE) [8,9], 
which accurately takes into account the transverse dy- 
namics for cigar-like condensates. We analyze the ex- 
istence, stability and collective oscillations of these 3D 
Bose-condensed bright solitons starting from a BEC with 
transverse confinement. We also study the existence and 
stability of multi-soliton condensates in a box and sim- 
ulate the dynamics of the " soliton train" experimentally 
observed in Ref. [5]. Then we investigate the macro- 
scopic quantum tunneling of a BEC on a Gaussian barrier 
showing that the bright-soliton strongly reduces interfer- 



ence fringes in the tunneling region. We show that the 
collapse of BEC can be induced by its scattering on a 
barrier if its density reaches a critical value, which can 
be analytically predicted. Finally we study the scatter- 
ing of bright solitons and the conditions for their collapse 
at the collision. 



II. BRIGHT SOLITONS UNDER TRANSVERSE 
CONFINEMENT 

The dynamics of a BEC at zero temperature is well de- 
scribed by the 3D Gross-Pitaevskii equation (3D GPE) 
[7] . In many cases the numerical solution of 3D GPE is a 
hard task and approximate procedures are needed. Start- 
ing from the 3D GPE we have recently derived and stud- 
ied [8,9] an effective time-dependent ID nonpolynomial 
nonlinear Schrodinger equation (NPSE). NPSE describes 
very accurately Bose condensates confined by a harmonic 
potential with frequency uj± and harmonic length a± = 
{Ti/moj^) 1 / 2 in the transverse direction and by a generic 
potential V(z) in the axial one. The total wave function 
of the condensate is ip( x i Vi z i t) = t)(j>(x, y, t), where 
the transverse wavefunction cf>(x,y,t) is a Gaussian with 
a width rj given by rj 2 = a\ y/l + 2a s iV|/| 2 and f(z,t) 



satisfies the NPSE equation 
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where a s is the s-wave scattering length. N is the 
number of condensed bosons and the function f(z,t) 
is normalized to one. In the weakly-interacting limit 
a s iV|/| 2 << 1, NPSE reduces to a ID GPE. Instead, 
in the strongly-interacting limit, NPSE becomes a non- 
linear Schodinger equation with the nonlinear term pro- 
portional to |/|/ [8,9]. 

Under transverse confinement and negative scattering 
length (a s < 0) a bright soliton sets up when the negative 
inter-atomic energy of the BEC compensates the positive 
kinetic energy such that the BEC is self-trapped in the 
axial direction. The shape of this 3D Bose-condensed 
bright soliton can be deduced from NPSE [8]. Setting 
V(z) = 0, scaling z in units of «i and t in units of wj , 
with the position 
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one finds the bright-soliton solution written in implicit 
form 
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FIG. 1. Axial density profile p(z) of the Bose-condensed 
bright soliton: 3D GPE (full line), Eq. 3 (dotted line), Eq. 5 
(dashed line). Length in units a± = (R/mc^x) 1 ^ 2 and density 
in units l/a±. 
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where ( = z — vt and 7 = \a s \N/a±. This equation is well 
defined only for 7$ 2 < 1/2; at 7<I> 2 = 1/2 the transverse 
size is zero. Moreover, by imposing the normalization 
condition one has 
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The normalization relates the chemical potential /1 to 
the coupling constant 7, while the velocity v of the 
bright soliton remains arbitrary. In the weak-coupling 
limit (7 ( E >2 << 1), the normalization condition gives 
fi = 1 — 7 2 /2 and the bright-soliton solution reads 
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The above solution is the text-book ID bright soliton 
of the ID nonlinear (cubic) Schrodinger equation (ID 
GPE). 

In Figure 1 the "exact" solitonic solution obtained by 
numerically solving the 3D GPE with a finite-difference 
predictor-corrector method [10] is compared with the an- 
alytical solution (Eq. 3) of the NPSE and the analytical 
solution (Eq. 5) of the ID bright soliton. The density 
profile is plotted for increasing values of 7. The agree- 
ment between the "exact" solution and the analytical one 
obtained from the NPSE is always remarkably good. In- 
stead, as expected, the profile of the ID bright soliton 
deviates from the 3D one for large 7. 



From Eq. 4 it easy to show that for 7 > 2/3 there 
are no solitary-wave solutions. Thus the NPSE gives the 
condition 
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for the existence of the 3D Bose-condensed bright soli- 
ton under transverse confinement. Note however that for 
7 C = 2/3 the transverse size of the Bose-condensed soli- 
ton is not zero, in fact 7 C $ 2 < 1/2. We have numerically 
verified by solving the 3D GPSE that the critical value 
7 C = 2/3 is very accurate. This is a remarkable result 
because, contrary to the 3D bright soliton (Eq. 3), the 
widely studied ID bright soliton (Eq. 5) exists (and it is 
stable) at any 7: for large values of 7 the wavefunction 
simply becomes narrower. 

In order to analyze the stability of our 3D soliton so- 
lution and its collective modes, it is useful to provide a 
simple analytical representation of its shape. The most 
natural choice is, as usual, a Gaussian with two varia- 
tional parameters: the transverse width a and the axial 
width 77 [11]. In this case, the energy per particle E of 
the condensed atomic-cloud obtained from the 3D GPE 
energy functional is given by 
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where g = y/2/irj, lengths are in units a± and the energy 
is in units Huj±. 
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FIG. 2. Gaussian approximation of the 3D bright soliton. 
Top: Widths a and tjasa function of the interaction strength 
7. Bottom: Collective frequencies wi and u>2 as a function of 
the interaction strength 7, with 7 = N\a s \/a±. Length in 
units a± = (S/mui) 1 ' 2 , time in units wj 1 . 

The energy function E = E(a, rj) has a local minimum, 
that is the 3D bright soliton condition, if the inter-atomic 
strength is smaller than a critical value: for larger values 
the soliton collapses. In Figure 2 we plot the widths a and 
r\ of the approximated 3D bright soliton as a function of 
the inter-atomic strength. Note that the transverse width 
a does not change very much while the longitudinal width 
7] is divergent for 7 = and it is comparable with a for 7 
close to 7 C . Moreover the critical value j c of the collapse 
is slightly overestimated with respect to the analytical 
prediction 7c = 2/3 of NPSE. 

The Gaussian approximation of the BEC wavcfunction 
can be used to study the dynamical stability of the 3D 
bright soliton by means of its collective oscillations. The 
diagonalization of the Hessian matrix of the energy func- 
tion E(a, rj) gives two frequencies lo\ and u>2 of collective 
excitations around the bright soliton solution. In Figure 
2 we plot such frequencies and the widths a and rj of 
the solitonic configuration as a function of the strength 
7. The frequencies are real and it means that the bright 
soliton is dynamically stable until it collapses at j c . Note 
that only for 7 = the frequencies ui\ and ui2 can be in- 
terpreted as transverse and axial collective oscillations of 
the bright soliton; nevertheless the mixing angle remains 
quite small also for finite values of 7 so they can be asso- 
ciated the transverse and axial motion respectively. The 
transverse frequency w\ is practically constant but close 
to 7 C it suddenly grows. The axial frequency L02 is zero 
for 7 = it increases with 7 but close to 7 C it goes to 
zero. 



III. MULTI-SOLITON BOSE CONDENSATES IN 
A BOX 

The solitary bright-soliton solution (Eq. 3) has been 
found by using Eq. 1 and Eq. 2. From these equa- 
tions one obtains the Newtonian second-order differential 
equation 
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where ( = z — vt and 7 = \a s \N/a±. The constant of 
motion of this equation is given by 



(9) 



By imposing the boundary condition $ — > for ( — > 00, 
which implies that E = 0, one has the solitary bright- 
soliton solution, that has only one peak. 




FIG. 3. Axial density profile p(Q of the bright 

multi-soliton Bose-condensate in a Box of length L = 10. 
From top to bottom: one, three and twelve solitons close to 
their critical stength. Units as in Fig. 1. 

It is important to observe that it is possible to obtain 
a multi-soliton static solution, i.e. a multi-peak solu- 
tion, by imposing box boundary conditions: L/2) = 
$(L/2) = 0, where L is the length of the axial box. The 
density profile is periodic and the quantized number N s 
of peaks defines the number of solitons in the box for 
a given chemical potential [i and inter-atomic strength 
7. In Figure 3 we plot the axial density profile obtained 
from the numerical integration of Eq. 9 for one, three 
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and twelve bright solitons in a box of length L = 10aj_. 
Also in the case of 3D bright multi-soliton Bose conden- 
sates, due to the finite transverse confinement, the solu- 
tion exists only if the interaction strength 7 = N\a s \/a± 
is below a critical value. From the normalization condi- 
tion it is easy to find that the condition of existence for 
a multi-soliton solution is given by 

L 



< Na s < 
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in the limit N s — > 00. Thus, in a large box, the strength 
7 has the critical value 7 C = 2/3 for one soliton and 
7 C = L/(4a±) for many-solitons. 
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Table 1. Bosc-Einstein condensate with N s bright soli- 
tons in a box of length L = 10a^. Chemical potential 
H and maximum density po of the multi-solitons at the 
critical strength j c . Units as in Fig. 1. 

In Table 1 we show the chemical potential fi, the max- 
imum density p and the critical strength j c of a Bose 
condensate made by N s bright solitons in a box of length 
L = 10a^. As expected the critical strength j c grows 
with N s and for large N s it approaches L/(4a±) — 2.5. In 
the case N s even, the multi-soliton solutions are also so- 
lutions with periodic boundary conditions (toroidal con- 
figuration with axial length L). In a torus, this train 
of solitons can travel without spreading with arbitrary 
velocity. 



IV. PROPAGATION OF A SOLITON TRAIN: 
COMPARISION WITH EXPERIMENTS 

In Rcf. [5] it has been reported the formation of a 
multi-soliton Bose condensate of 7 Li atoms created in 



a quasi- ID optical trap that can be modelled by an 
asymmetric harmonic potential with a large aspect ra- 
tio (A = u)±/uj z = 100). In this experiment, a "soli- 
ton train", containing N s solitons (N s = 4,5,6), has 
been created from a Bosc-Einstein condensate by mag- 
netically tuning the inter-atomic interaction from repul- 
sive (a s = Was, with as the Bohr radius) to attractive 
(a s = — 3as) using a Fcshbach resonance. The soliton 
train has been set in motion and observed to propagate 
in the axial optical harmonic potential for many oscilla- 
tory cycles with a small spreading. Moreover, it has been 
found that the spacing between solitons is compressed 
at the turning points and spread out at the center of 
the oscillation, suggesting a repulsive interaction between 
neighboring solitons. 
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FIG. 4. Axial density profile of the travelling "soliton 
train" of 7 Li atoms in a anysotropic harmonic potential with 
uj± — 400 Hz and lu z = 4 Hz. Scattering length: a s = San 
with as Bohr radius. 5000 atoms for each soliton. Units as 
in Fig. 1 and Fig. 2. 

By using our time-dependent NPSE it is quite easy to 
simulate the dynamics of the soliton train described in 
[5]. In Figure 4 we plot three frames of the time evo- 
lution of a set of four 7 Li bright solitons, each of them 
made of 5000 atoms, initially located far from the cen- 
ter of the axial harmonic potential with an alternating 
phase structure. Note that the train configuration ceases 
to be a true solitonic solution as soon as we switch off the 
axial confinement. In fact its shape and density change 
during the motion. The agreement between our numeri- 
cal results of Figure 4 and the experimental data shown 
in Figure 4 of Ref. [5] is remarkable: the spacing be- 
tween solitons increases near the center of oscillation and 
bunches at the end points. An alternating phase struc- 
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ture is essential to have a repulsion between neighboring 
solitons [12]. In fact, we have verified that a train of soli- 
tons with the same phase (attractive interaction between 
solitions) becomes a narrow blob of high density which 
can induce the collapse of the condensate. In the case 
of positive scattering length our calculations show that 
during the motion an initially confined Bosc condensate 
made of many peaks becomes a blob which spreads along 
the axial trap but partially recomposes at the turning 
points. 

The effective interaction between two bosonic matter 
waves does indeed depend on their phase difference A , 
being attractive for A > and repulsive otherwise. This 
phenomenon is fully analogous to what has been already 
studied in the framework of non linear optics and may 
be simply understood on the basis of an elementary argu- 
ment regarding the extent of the overlap between the two 
waves. From the one dimensional GPE, the interaction 
energy density of the condensate is proportional to the 
square of the local particle density: W(z) oc a s p(z) 2 . 
Now we consider a wave function accurately approxi- 
mated by the superposition of two well separated solitons 
%l>{z) ~ [<f>(z - z )e iA / 2 + <j>(z + z )e- lA ' 2 ] /y/2, where 
2zo is the relative separation. By substituting this form 
into the expression of the interaction energy density at 
z = we get H^(0) oc a s (cosA + 1). A measure of 
the effective interaction can be obtained by considering 
the difference between W^(0) and the analogous self in- 
teraction density of a single soliton Wi(0) which gives: 
W 2 (0) — Wi(0) <x a s cos A. In the case of negative scat- 
tering length the interaction turns out to be attractive 
as soon as cos A > and repulsive otherwise, in agree- 
ment with more quantitative analysis [12] and previous 
numerical findings. 

V. TUNNELING WITH BRIGHT SOLITONS 

Now we investigate the role of solitonic configurations 
and the effect of an attractive inter-atomic interaction on 
the behavior of a Bose-Einstein condensate during tun- 
neling. The initial condition of the BEC is the ground- 
state of NPSE with a harmonic trapping potential also in 
the horizontal axial direction: V(z) = ^muj 2 (z — zq) 2 . To 
have a cigar-shaped condensate we choose A = ljj_/u> z — 
10. We set z = 20, where z is written in units of the 
harmonic length a z = {Ti/muiz) 1 / 2 . For t > the trap in 
the axial direction is switched off and a Gaussian energy 
barrier is inserted at z = 0. The potential barrier is given 

by 

V(z) = V e~^' , (11) 

where Vq is the height of the potential barrier and S its 
width. The BEC is moved towards the barrier by adding 
an initial momentum pq in the axial direction: 
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FIG. 5. Axial density profile of the Bose condensate tun- 
neling through the Gaussian barrier. Start-up kinetic energy 
per particle of the condensate: Eq — p§/(2m) = 10. Gaus- 
sian barrier parameters: Vo = 10 and E = 1. Three values 
of the interaction strength: Na s /a z — —2 ■ 10 _1 (solid line), 
Na s /a z = (dashed line), Na s /a z = 2 ■ 10 _1 (dotted line). 
Length in units a z = (h/mtUz) 1 ^ 2 , time in units wj 1 , and 
energy in units hui z . 

We have verified that the condensate with a scatter- 
ing length close to the collapse value (Na s /a z = —0.2) 
is a bright soliton: its shape does not depend on the 
axial harmonic trapping potential. In Figure 5 we plot 
the axial density profile of the Bose condensate at differ- 
ent instants for three values of the interaction strength: 
Na s /a z = 0.2, Na s /az = and Na s /a z = -0.2. While 
the case of positive scattering length is quite similar to 
that of zero scattering length, the case with negative scat- 
tering length is peculiar: as previously stated the conden- 
sate is a soliton, moreover it does not show interference 
patterns and simply splits into a reflected and a transmit- 
ted matter wave. However, the reflected and the trans- 
mitted waves are not solitonic, i.e. the tunneling through 
a barrier breaks the solitonic nature of the matter wave. 

If we decrease further the scattering length (a s < 0), 
the Bose condensate becomes unstable. For Na s /a z — 
—0.2 the Bose-condensed bright soliton is still stable but 
we can induce its collapse at the impact time by increas- 
ing the initial momentum po or the energy barrier Vq . In 
fact, in these cases the local density of the condensate 
becomes sufficiently large to give rise to the collapse of 
the condensate. An inspection of Eq. 1 shows that a 
sufficient condition for the collapse of the condensate is 
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that 



«^l/l 2 < --. 



(13) 



that is the condition for which condensate shrinks to 
zero in the transverse direction. Our numerical compu- 
tations show that this condition gives precisely the one- 
dimensional density at which there is the collapse of the 
condensate during tunneling: pm = l/(2|a s |). 



VI. COLLISIONS WITH BRIGHT SOLITONS 

Solitons show peculiar properties not only in the col- 
lision with a barrier but also in the collision with other 
waves. In Figure 6 we plot some frames of the collision of 
two Bose-condensed bright solitons obtained by numeri- 
cally solving the time-dependent NPSE (Eq. 1). 
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FIG. 6. Collision of two bright solitions with a different 
number of particles: Ni = 2N2 and N\\a 8 \/a z = 0.1. Start-up 
momentum: po = 50. Units as in Fig. 5. 

The two solitons have a different number of particles 
and start with opposite momenta of equal modulus p . 
As shown in Figure 6, at the impact many fringes of inter- 
ference are produced. We have verified that the number 
of fringes grows with the momentum po while the max- 
imal density of the interference peak remains constant. 
Moreover, we have found that the only effect of a phase- 
difference between the two colliding solitons is a shift in 
the position of fringes. In particular, at a fixed time, by 
changing the phase difference by tt the spatial locations 
of maxima of interference are shifted into the positions 
of minima. After the impact the two bright solitons con- 
tinue their motion without any recollection of the im- 
pact: the solitons acquire again their initial shape and 
then continue the motion without shape deformations. 

This phenomenon of transparency can be also seen 
in the collision between a Bose-condensed bright soliton 



and a non-solitonic Bose condensate while two colliding 
generic matter- waves are not transparent. The numeri- 
cal solution of NPSE confirms that, apart the interfer- 
ence at the impact, the solitonic matter- wave and the 
non-solitonic matter-wave do not see each other. At large 
times the bright soliton moves with its initial density pro- 
file while the non-solitonic wave evolves with a spreading 
that is not influenced by the collision with the bright 
soliton. 

In the collision process of two Bose-condensed bright 
solitons the maximal density of the interference peak de- 
pends on the interaction strength, i.e. on the initial den- 
sity of the two colliding solitons. For two colliding bright 
solitons with the same number of particles the maximal 
density at impact is about four times the initial density 
of each soliton. We have verified that when the maximal 
density satisfies the condition of Eq. 13 then the collapse 
of the condensates at the impact sets in. 



CONCLUSIONS 

Our calculations suggest that, contrary to the ID 
bright soliton, the 3D bright soliton under transverse 
confinement exists and it is dynamically stable only if 
the attractive inter-atomic interaction is smaller than a 
critical value. We have analytically determined this crit- 
ical interaction strength and the density profile of the 3D 
bright soliton. We have also found the collective oscilla- 
tions of the bright soliton as a function of the interac- 
tion strength by using a Gaussian approximation. Then 
we have investigated multi-soliton Bose condensates with 
box boundary conditions. Also these multi-soliton solu- 
tions exist if the interaction strength is below a critical 
threshold that grows with box length. The dynamical 
properties of a train of solitons have been investigated 
by numerically solving our effective ID nonpolynomial 
nonlinear Schrodingcr equation: the results are in quali- 
tative agreement with the experimental data. In the case 
of soliton scattering, our theoretical results show that the 
collapse of 3D soliton can be induced by the scattering 
on a Gaussian barrier when the density at the impact 
reaches a critical value, which does not depend on the 
solitonic nature of the incident wave. A clear signature 
of solitonic behavior is the transparency during collisions. 
Analyzing the collision between bright solitons we have 
verified that, apart the interference at the impact, the 
Bose-condensed soliton is transparent, namely at large 
times it recovers its initial density profile and then travels 
without spreading. However, if the density of the interfer- 
ence peak exceeds a critical density the system collapses. 
Moveover we have shown that there is not transparency 
during tunneling: a bright soliton tunneling through a 
barrier splits into reflected and transmitted wave packets 
which are not solitonic. 



G 



[I] M.H. Anderson et al, Science 269, 189 (1995); K.B. 

Davis et al, Phys. Rev. Lett. 75, 3969 (1995); 
C.C. Bradley et al, Phys. Rev. Lett. 75, 1687 
(1995). 

[2] Bose-Einstein Condensation in Atomic Gases, vol. 
140 International School of Physics Enrico Fermi, 
Eds. M. Inguscio, S. Stringari, C. Wieman (IOS 
Press, Amsterdam, 1999). 

[3] P.G. Drazin and R.S. Johnson, Solitons: An Intro- 
duction (Cambridge University Press, Cambridge, 
1988). 

[4] S. Burger et al., Phys. Rev. lett. 83, 5198 (1999); 
J. Denshlag et al, Science 287 97 (2000). 

[5] K.E. Strecker et al, Nature 417, 150 (2002). 

[6] W.P. Reinhardt and C.W. Clark, J. Phys. B 30, 
L785 (1997); V.M. Perez-Garcia, H. Michinel, and 
H. Herrero, Phys. Rev. A 57, 3837 (1998); Th. 
Busch and J. R. Anglin, Phys. Rev. Lett. 87 
010401 (2001). 

[7] E.P. Gross, Nuovo Cimcnto 20, 454 (1961); L.P. 
Pitacvskii, Zh. Eksp. Teor. Fiz. 40, 646 (1961) 
[English Transl. Sov. Phys. JETP 13, 451 (1961)]. 

[8] L. Salasnich, A. Parola, and L. Reatto, Phys. Rev. 
A 65, 043614 (2002). 

[9] L. Salasnich, Laser Physics 14, 198 (2002). 

[10] L. Salasnich, A. Parola, and L. Reatto, Phys. Rev. 
A 64, 023601 (2001). 

[II] L. Salasnich, Int. J. Mod. Phys. B 14, 1 (2000). 
[12] J.P. Gordon, Opt. Lett. 8, 596 (1983). 



7 



